DETERMINAREA SOLUŢIILOR ECUAŢIILOR DIOFANTICE 

2069-2080 

CIRA OCTAVIAN 



REZUMAT. în acest articol ne propunem să determinăm soluţii 
pentru ecuaţiile diof antice 2069-2080. 



1. Introducere 

Funcţia care asociază fiecărui număr natural n pe cel mai mic nu- 
măr natural m care are proprietatea că ml este multiplu lui n a fost 
considerată de prima dată de Lucas, (3), în anul 1883. Alţi autori care 
au considerat această funcţie în lucrările lor sunt: Nueberg (1887) 0], 
Kempner (1918) |2|- Această funcţie a fost redescoperită de Smaran- 
dache, 85J, în anul 1980. Funcţia este notată de Smarandache cu S 
sau cu 77 , iar în site-ul Wolfram Math World, [ŢT 0 |, se notează cu fi. In 
acest articol am adoptat notaţia r; a lui Smarandache din lucrarea [[ZD- 
In monografia UZI, Smarandache prezintă o listă de ecuaţii diofan- 
tice referitoare la funcţia r/. Fie m, n e N* fixaţi şi x. y e N necunos- 
cutele. Ecuaţiile diofantice 2069-2080 cu funcţia r] sunt: 

(2069) 77(777 ■ x + n ) = x , 

(2070) 77(777 ■ x + n) — m + n ■ x , 

(2071) rj(m ■ x + n) = x\ , 

(2072) r]\x m ) = x n , 

(2073) r){x) m = r](x n ) , 

(2074) rj(rn ■ x + n) = r](x) y , 

(2075) Tj(x) + y = x + r)(y) , unde x şi y nu sunt prime. 

(2076) rj(x) + îj{y) = y{x + y), unde x şi y nu sunt prime gemene. In 
general funcţia 77 nu este o funcţie aditivă. 

(2077) rj(x + y) = rj(x) ■ r](y) . în general funcţia 77 nu este o funcţie 
exponenţială. 

(2078) r](x ■ y) — 77 (x) • 77 ( 7 /) . în general funcţia 77 nu este o funcţie 
multiplicativă. 

(2079) 77(771 • x + n ) = x y , 

(2080) rj(x) ■ y = x ■ 77 ( 7 /), unde x şi y nu sunt prime. 
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2. Algoritmul Kumpner 

Funcţia rj : N* — » N*, 77(71) = m, unde m este cel mai mic număr 
natural care are proprietatea că n îl divide pe m\, (sau m\ este mul- 
tiplu de n) este cunoscută în literatura de specialitate ca funcţia lui 
Smarandache. Valorile funcţiei, pentru n — 1 , 2, . . . , 18 , sunt: 1 , 2, 3 , 4 , 
5 , 3 , 7 , 4 , 6, 5 , 11 , 4 , 13 , 7 , 5 , 6, 17 , 6 obţinute cu un algoritm ce rezultă 
din definiţia funcţiei //, după cum urmează: 



Programul 2.1. 



di 71 ) 



for m = l..n 

return m if mod(m\,n ) = 0 



Acest algoritm nu poate fi folosit pentru n > 19 pentru că nume- 
rele 19 !, 20 !, . . . sunt numere cu mai mult de 17 cifre zecimale şi în 
modul clasic de calcul (fără o aritmetică a preciziilor arbitrare 0) se 
vor genera erori de trunchiere. 

Kempner, O, a dat un algoritm pentru a calcula pe 77(71) folosind 
factorizarea clasică pff ■ pff ■ ■ ■ pf s cu numere prime a lui n, şi baza 
de numeraţie generalizată (n,;)^], pentru i = l,s. Soluţii parţiale 
pentru algoritmul de calcul a lui 77(71) s-au dat mai devreme de Lucas 
şi Neuberg, 0 • 

Prezentăm algoritmul lui Kempner, J 2 J, de calcul al funcţiei lui 
Smarandache 77. Pentru început definim şirul recursiv a J+i = p ■ a 3 + 1 
cu j = 1,2 ,... şi ai = 1 , unde p este un număr prim. Acest şir de 
numere constituie baza de numeraţie generalizată a lui p. Deoarece 
a 2 = p + 1 , a 3 = p 2 + p + 1 , ... se poate demonstra prin inducţie că 
dj — l+p + ... + p j ~ 1 = {jf — l)/{p — 1 ) pentru Vj > 2. 

Valoarea lui v, astfel încât a v < a < a u+ 1, este dat de formula 



(2.1) v = [logp (1 + a(p - 1))J , 

unde L'J este funcţia parte întreagă inferioară. Cu ajutorul algoritmului 
lui Euclid putem determina şirurile unice Kt şi r, după cum urmează 



(2.2) 


a = 


^iz * ~\ ~ V* 1/ 1 


(2.3) 


A = 


k v -\ • 1 + r^_i , 


(2.4) 






(2.5) 


T v-{ A-2) = 


(A— 1) * &v— (A— 1) “ 1 “ ^ v— (A— 1) 


(2.6) 


A-(a-i) = 


1 

V 

P 

1 
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adică, până când restul r„_ \ = 0. La fiecare pas este partea întreagă 
a raportului rr/a, şi r, este restul împărţirii. De exemplu pentru pri- 
mul pas avem k v = [a/a v \ şi r u = a — k v ■ a v . Atunci avem 

A 

(2.7) rj(p a ) = (p - l)a + ^ . 

i=v 

în general pentru 

( 2 . 8 ) n = p? -pT-'-pT, 

valoarea funcţie r; este dată de formula: 

(2.9) p(n) = max {^(p? 1 ), v(p?), • • ■ , v{pT)} , 
formulă ce o datorăm lui Kempner |2|. 



Remarca 2.2. Dacă n e N* are descompunerea în produs de numere 
prime (2.8), unde pi sunt numere prime astfel încât pi < p 2 <■■■ < 
p s , iar s > 1, atunci algoritmul lui Kempner de calcul al funcţiei // 
este 



(2.10) V {n) = max |p, • (oi [ pi] ) , p 2 ■ (% 2] ) fc) , 

Ps ' ’ 

unde prin (<%]) înţelegem că numărul a este "scris" în baza de nu- 
meraţie generalizată p (notată aq p ]) şi este "citit" în baza de numeraţie 
p (notat /3( p ), unde Ş = U[ p ]), B p. 39]. 

3. Programe 



Programul 3.1. Funcţia care numără cifrele în baza de numeraţie p a 
lui n 



ncb(n,p ) 



return ceil(\og(n,p)) if n > 1 
return 1 otherwise 



unde funcţia utilitară Mathcad ceil(-) este funcţia parte întreagă su- 
perioară. 



Programul 3.2. Program pentru generarea a bazei de numeraţie ge- 
neralizată a lui p (notată de Smarandache [p]) pentru un număr cu 
m cifre 



a(p, m) 



for i E l..m 

p l — 1 

di ^ Ţ- 

p — 1 

return a 
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Programul 3.3. Program de generarea a bazei de numeraţie p (notată 
de Smarandache (p)) pentru a putea scrie numărul a 



b(at,p ) 



return (1) if p = 1 
for i G 1 ..ncb(a,p) 
bi G- p * _1 
return b 



Programul 3.4. Programul determină cifrele bazei de numeraţie ge- 
neralizată [p] pentru numărul n 



Nbg(n,p) := 



m G- ncb(n,p) 
a G- a(p, m) 
return (1) if m — 0 
for 



n 



c,; G- trunc | — 

di 

n G- mod (n, a*) 
return c 

Programul 3.5. Program pentru funcţia Smarandache 

77(71) := return ” Err. n nu este întreg'’’’ if n ^ trunc(n ) 
return ” Err. n < 1 ” if n < 1 
return (1) if n = 1 
/ G- Fain) 

P / (1> 
o G- / <2> 

/or k = l..rows(p ) 

Vk G- p k ■ Nbg(a kj p k ) • b(a k ,Pk ) 
return max(p) 

Acest program apelează programul Fain) de factorizare cu numere 
prime. Programul foloseşte remarca 2.2 a lui Smarandache referitor 
la algoritmul Kempner. 

Programul 3.6. Programul de generarea numerelor prime până la L 
cu ajutorul ciurului lui Atkin. 



Atkin(L) : = 



eprirriL G- 0 
A G- floor(g r L) 
for j G 1..A 



for k G 1. .cei/ 



v/r 






;2 



71 G- 4fc 2 + J 2 
777 G- r mod(n, 12) 

eprim n G- ~^eprim n ifn<LA(m — l\/m — 5) 
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for k G 1. .ceil 






n •<— 3 k 2 + j 2 

eprim n •<— ->eprim n if n < L A mod(n, 12) = 7 
/or k e j + \..ceil ' A+T 
n •<— 3/c 2 — j 2 

eprim n •<— -<eprim n if n < L A mod(n, 12) = 11 
for k G 5, 7, ..A 

forj — 1, 3, if eprim k 

eprirrij . k 2 0 
eprirrii <— 2 
eprim 2 G- 3 
j*— 3 

/or n G 3,5,..L 
if eprim n 
prime j •<— n 
J A j + 1 
return prime 



Programul 3.7. Programul de factorizare a unui număr natural. 

Acest program foloseşte vectorul numerelor prime p care este ge- 
nerat cu ciurul lui Atkin, na, cel mai rapid program de generare a 
numerelor prime până la limita precizată. Apelul ciurului lui Atkin 
se face cu secvenţa: 

fo = time{ 0) 



p := Af/cm(10 8 ) Zasf(p) = 5761455 p Zast(p) = 99999989 



fi = fime(l) (fi — t 0 )sec = 82.88 sec . 



Fa(m) 



return (”m = ” m ” > ca ultimul p 2n ) 

j 1 

k <- 0 



while m > p.,- 

if mod ( m,pj ) = 0 

k <— k + 1 
m 

mi 

Pj 

otherwise 



if m > (piaspp)) 2 
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/ A- stack[f, ( pj , /c)] if k > 0 
J «- J + 1 
/c A- 0 

/ A- stac/c[/, (pj, k)\ if k > 0 
return submatrix(f , 2, rows(f), 1 , 2) 

Programul 3.8. Program de calcul al valorilor funcţiei r/ pentru nume- 
rele naturale până la IO 6 . 

ValFS(N) := VFS 1 A- 1 
for n E 2..N 
' VFS n A- 77 ( 71 ) 
return V FS 

Apelul acestui program va genera un vector care va fi salvat ca un 
fişier de tip prn. Astfel vom putea folosi acest fişier pentru rezolvarea 
ecuaţiilor diofantice. Cu ajutorul funcţie Mathcad tirrie s-a contorizat 
timpul de generare al vectorului VFS. 

N := IO 6 t 0 = time{ 0) VFS := ValFS(N) t x : time{ 1) 

(ti — t 0 ) ■ sec = ”1:7: 32.623” • hhmmss 
W RIT EPRN (” V FS.prn” ) := VFS last(VFS ) = 1000000 
Pentru informare se afişează ultimul indice al vectorului VFS. 

4. Rezolvarea ecuaţiilor diofantice 

4.1. Ecuaţia (2069). Pentru căutarea soluţiilor vom citi fişierul V FS.prn 
cu ajutorul funcţiei Mathcad READPRN 

VFS := RE ADPRN(” ...\V FS.prn ”) last(VFS ) = 1000000 

unde comanda last(V SF) ne indică ultimul indice al vectorului V FS. 

Generarea fişierului V FS odată şi citirea fişierului generat în do- 
cumente Mathcad care determină soluţii ale ecuaţiilor diofantice face 
ca să rezulte o mare economie de timp de execuţie a programului de 
căutare a soluţiilor. 

Programul 4.1. Programul de determinare a soluţiilor ecuaţiei (2069). 

Ed2069(a m , b m , a n , b n , a x , b x ) := S A- (”m” "n" ”x") 

u E- last(VFS) 
for m G a m ..b m 
for ne a n ..b n 
for x e a x ..b x 
rj A- m ■ x + n 
q A < u A V FS^ = x 
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S G- stack[S, (m n x)] if q 
return S 

Apelul programului se face prin secvenţa: 

ci m . 2 b m . 10 a n . 1 b n . 10 a x . 1 b x . 16 

Numărul total de cazuri verificate este: 

{bm — + 1) {bn ~ 0>n + 1 ){f>x ~ a x + 1) = 1440 . 

Apelul programului Ed2069: 

t 0 : time{ 0) Sol := Ed2069{a m , b m , a n , b n , a x , b x ) ti := time{ 1) 

Timul de execuţie în secunde şi numărul de soluţii rezultă din: 

{ti — to) ■ s = 0.011 • s rows{Sol ) — 1 = 36 

Pentru m G (2, 3, , 10}, n G (1, 2, ... , 10} şi x G (1,2,..., 16} cele 
36 de soluţii ale ecuaţiei diofantice rj{m • x + n) = x sunt: 




Valoarea maximă a soluţiilor x este 15. 

4.2. Ecuaţia (2070). Pentru căutarea soluţiilor vom citi fişierul V FS.prn 
cu ajutorul funcţiei Mathcad READPRN 

VFS := READPRN (” . . \V FS.prn ” ) last{VFS) = 1000000 . 

Programul 4.2. Programul de determinare a soluţiilor ecuaţiei (2070). 

Ed2070{a m , b m , a n , b n , a x , b x ) := S G- (”m” ”n” ”x”) 

u •(— last{VFS) 
for m G a m ..b m 
for n G a n ..b n 
for x G a x ..b x 
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rj E- m ■ x + n 

q <— V < u A V FS v — m + n ■ x 
S E- stack[S, (m n x)] if q 
return S 

Apelul programului se face prin secvenţa: 

a m := 2 b rn := 20 a n := 1 b n := 20 a x := 1 b x := 16 
Numărul total de cazuri verificate este: 

(b m — o>m + 1 )(b n — a n + 1 ) ipx — a x + 1) = 7220 . 

Apelul programului Ed2070: 

to : time( 0) Sol := Ed2070(a m , b m , a n , b n , a x , b x ) t\ := time( 1) 

Timul de execuţie în secunde şi numărul de soluţii rezultă din: 

(ti ~ to) ■ s = 0.853 • s rows(Sol) — 1 = 14 

Pentru m G (2, 3, , 20}, n E (1, 2, ... , 20} şi x E (2, 3, ... , 20} cele 
14 de soluţii ale ecuaţiei diofantice r/(m • x + n) = x sunt: 



m 


2 


4 


4 


6 


6 


8 


10 


12 


12 


14 


16 


18 


18 


20 


n 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


X 


4 


2 


6 


4 


8 


6 


12 


10 


14 


12 


18 


16 


20 


18 



Valoarea maximă a soluţiilor x este 20. 

4.3. Ecuaţia (2071). Pentru căutarea soluţiilor vom citi fişierul V FS.prn 
cu ajutorul funcţiei Mathcad READPRN 

VFS := RE ADPRN(”...\V FS.prn”) lastiVFS) = 1000000 . 

Programul 4.3. Programul de determinare a soluţiilor ecuaţiei (2071). 

Ed207l(a m ,b m ,a n ,b n ,a x ,b x ) S <— (”m” ”n” ”x”) 

u E- last(VFS) 
for m E a m ..b m 
for n E a n ..b n 
for x E a x ..b x 
rj E- m ■ x + n 
q E- ij < u A VFS^ = x\ 

S E- stack[S, (m n x)] if q 
return S 

Apelul programului se face prin secvenţa: 

a m . 2 b m . 15 a n . 1 b n . 15 a x . 1 b x . 19 
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Numărul total de cazuri verificate este: 

( b m — a m + 1 )(b n — a n + 1 )(b x — a x + 1) = 3990 . 

Apelul programului Ed2071: 

to : time( 0) Sol := Ed2071(a m , b m , a n , b n , a x , b x ) t\ := tvme( 1) 

Timul de execuţie în secunde şi numărul de soluţii rezultă din: 

(ti — to) ■ s — 0.02 • s rows(Sol) — 1 = 24 

Pentru m G (2, 3, , 15}, n G (1, 2, ... , 15} şi x G (1,2,..., 19} cele 
24 de soluţii ale ecuaţiei diof antice 77(771 ■ x + n) = x sunt: 



/ 777 N 




f 2 


2 


2 


3 3 


4 


4 


5 5 


7 


8 


9 


10 




n 


= 


3 


10 


12 


7 9 


4 


6 


1 3 


15 


12 


9 


6 




V x > 






3 


3 


3 3 


3 


3 


3 3 


3 


3 


3 


3 








10 


11 


11 


11 


12 


12 


13 


13 


14 


14 


15 






15 


3 


12 


15 


9 


12 


6 


9 


3 


6 


3 






3 


3 


3 


3 


3 


3 


3 


3 


3 


3 


3 



Valoarea maximă a soluţiilor x este 3. 



4.4. Ecuaţia (2072). Pentru căutarea soluţiilor vom citi fişierul V FS.prn 
cu ajutorul funcţiei Mathcad READPRN 

VFS := RE ADPRN(”...\V FS.prn”) last(VFS) = 1000000 . 



Programul 4.4. Programul de determinare a soluţiilor ecuaţiei (2072). 

Ed2072(a m , b m , b n , (ix, b x ) ■ 



’m” ”77” ”x” 



u last(VFS) 
for 777 G a m ..b m 
for n G a n ..b n 
for x G a x ..b x 
77 G- 

g G- 77 < 77 A V = T n 
S 1 stac/cfS 1 , (777 77 t )] if q 
return S 



Apelul programului se face prin secvenţa: 

a m . 2 b m . 9 <7 n . 2 b n . 9 a x . 2 b x . 10 

Numărul total de cazuri verificate este: 

(bm ®rn 1 )(bn Ţ- 1 )(b x o x G 1) — 0/6 . 

Apelul programului Ed2072: 

to : time(0 ) Sol := Ed2072(a m , b m , a n , b n , a x , b x ) ti := time( 1) 
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Timul de execuţie în secunde şi numărul de soluţii rezultă din: 

(t\ — to) ■ s — 0.11 • s rows(Sol) — 1 = 12 

Pentru m G (2,3,..., 9}, n G (2, 3, . . . , 9} şi x G (1,2,..., 10} cele 12 
de soluţii ale ecuaţiei diofantice r/(x m ) = x n sunt: 



/ m \ / 2 


3 


3 


4 


5 


5 


6 


6 


6 


7 


7 


7 


| n | = |2 


2 


2 


2 


2 


3 


2 


2 


3 


2 


2 


3 


\ x ) \ 2 


2 


3 


3 


5 


2 


4 


5 


2 


4 


7 


2 



Valoarea maximă a soluţiilor x este 7. 

4.5. Ecuaţia (2073). Pentru căutarea soluţiilor vom citi fişierul V FS.prn 
cu ajutorul funcţiei Mathcad READPRN 



VFS := READPRN (” . . .\V FS.prn” ) last(VFS ) = 1000000 . 



Programul 4.5. Programul de determinare a soluţiilor ecuaţiei (2073). 



Ed2071( ®i») b m , b n , bx) 



(”m” ”n” ”x”) 
u •(— last(VFS ) 
for m G a m ..b m 
for n e a n ..b n 
for x G a x ..b x 

q^x n <uA ( VFS x ) m = VFS X , 
S A- stack[S, (m n a;)] if q 
return S 



Apelul programului se face prin secvenţa: 

a m . 2 b m . 9 a n . 2 b n . 9 o x . 2 b x . 25 

Numărul total de cazuri verificate este: 

(bm + l)(b n o n + 1 ~}(bx o x Al) 1536 . 

Apelul programului Ed2073: 

t 0 : time( 0) Sol := Ed2073(a m , b m , a n , b n , a x , b x ) ti := time( 1) 

Timul de execuţie în secunde şi numărul de soluţii rezultă din: 

(ti — t 0 ) ■ s = 0.014 • s rows(Sol) — 1 = 20 
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Pentru m e {2,3,..., 9}, n G {2, 3, . . . , 9} şi x G {2, 3, ... , 25} cele 20 
de soluţii ale ecuaţiei diofantice r/(x) m = rj(x n ) sunt: 



/ rn \ / 2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


| n J = I 2 


3 


3 


3 


4 


4 


4 


4 


5 


5 


5 


5 


6 


6 


\ x ) \ 2 


2 


3 


6 


3 


6 


8 


24 


5 


8 


10 


15 


4 


5 



2 2 2 3 3 3 \ 

6 7 7 5 6 7 

10 4 7 2 2 2 ) 

Valoarea maximă a soluţiilor x este 24. 
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